INSTITUT FUR KERNPHYSIK, UNIVERSITAT FRANKFURT 

D - 60486 Frankfurt, August-Euler-Strasse 6, Germany 



IKF-HENPG/2-99 



a^ ■ Evidence for Statistical Production of J/ip Mesons in 

Nuclear Collisions at the CERN SPS 

o 

oo 

^ . Marek GazdzickiQ 

! Institut fiir Kernphysik, Universitat Frankfurt, Germany 

in 

in \ Mark I. Gorenstein[|'Q 

O ■ Institut fiir Theoretische Physik, Universitat Frankfurt, Germany 

D , Abstract 

The hypothesis of statistical production of J/'0 mesons at hadronization is formu- 
lated and checked against experimental data. It explains in the natural way the 
r\ \ observed scaling behavior of the J/ijj to pion ratio at the CERN SPS energies. Us- 

' ing the multiplicities of J /ip and r] mesons the hadronization temperature Th = 170 

MeV is found, which agrees well with the previous estimates of the temperature 
parameter based on the analysis of the hadron yield systematics. 
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Charmonium production in hadronic |^ and nuclear |0] collisions is usually considered 
to be composed of three stages: the creation of a cc pair, the formation of a bound cc state 
and the subsequent interaction of this cc bound state with the surrounding matter. The 
first process is calculated within perturbative QCD, whereas modeling of non-perturbative 
dynamics is needed to describe the last two stages (see, e.g. and references therein). 
The interaction of the bound cc state with matter causes suppression of the finally ob- 
served charmonium yield relative to the initially created number of bound cc states. This 
initial number is assumed to be proportional to the number of Drell-Yan pairs, which 
then allows for the experimental study of the charmonium suppression pattern. It was 
proposed 0, |^ that the magnitude of the measured suppression in nuclear collisions can 
be used as a probe of the state of high density matter created at the early stage of the 
collision. The suppression of the J /ip yield observed in p+A and 0(S)+A collisions at the 
CERN SPS is considered to be caused by the interactions with nucleons occurring while 
the primordial baryons keep interpenetrating |]^. The rapid increase of the suppression 
{anomalous suppression) observed when going from peripheral to central Pb+Pb collisions 
is often attributed to the formation of a quark-gluon plasma However alternative 
interpretations are still under discussion p, |], ||. 

It was recently found [0, |ll| that the mean multiplicity of J / ip mesons increases pro- 
portionally to the mean multiplicity of pions when proton-proton (p+p), proton-nucleus 
(p+A) and nucleus-nucleus (A+A) collisions at CERN SPS energies are considered. We 
illustrate this unexpected experimental fact by reproducing in Fig. |l| the plot from Ref. 
||1 1|| , where the ratio {J/tp)/{h~) is shown as a function of the mean number of nucleons 
participating in the interaction for inelastic nuclear collisions at the CERN SPS. The 
(J/ip) and {h~) denote here the mean multiplicites of J/ip mesons and negatively charged 
hadrons (more than 90% are n~ mesons), respectively. We note that the analysis pre- 
sented in Ref. |]TD[ indicates that the scaling of the {J/ip)/{h~) ratio is also valid for 



central Pb+Pb collisions at the CERN SPS. 

In the standard picture of the J/ ip production based on the hard creation of cc pairs and 
the subsequent suppression of the bound cc states the observed scaling behavior of the 
J /ip multiplicity appears to be due to an 'accidental' cancelation of several large effects. 
This motivates our effort to find an alternative production mechanism of J /if) mesons 
which would explain the experimental data in a natural way. 
In this letter we show that a scaling property of the J/ip multiplicity 

^ const{A) (1) 

can be understood assuming that a dominant fraction of J/ip mesons is produced directly 
at hadronization according to the available hadronic phase space. 

Since a long time |T2[ statistical models are used to describe hadron multiplicities in 



high energy collisions. Thermal hadron production models have been successfully used 
to fit the data on particle multiplicities in A+A collisions at the CERN SPS energies 



(see, e.g. [13, |T^). Due to the large number of particles a grand canonical formulation is 
used for the modeling of high energy heavy ion collisions [jl5|, 0] . Recently, an impressive 
success of the statistical model applied to hadron multiplicities in elementary e"*" + e~. 
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p + p and p + p interactions at high energy was also reported |T^. However, in the latter 
case the use of a canonical formulation of the model, which assures exact conservation 
of the conserved charges, is necessary. The temperature parameter which characterizes 
the available phase space for the hadron production is found in these interactions to be 
160-190 MeV It does not show any significant dependence on the type of reaction 
and on the collision energy. Moreover, it coincides with the chemical freeze-out tem- 
perature estimate obtained in hadron gas models for A+A collisions at the CERN SPS 
||T3| . These facts suggest the possibility to ascribe the observed statistical properties of 
hadron production systematics in elementary and nuclear collisions at high energies to 



the statistical nature of the hadronization process ||16| , \L% [17 . 

Based on the above facts we formulate a hypothesis that a dominant fraction of 
the J/ip mesons produced in hadronic and nuclear collisions at the CERN SPS 
energies is created at hadronization according to the available hadronic phase 
space. 

J/ip mesons are neutral and unflavored, i.e. all charges conserved in the strong interac- 
tion (electric charge, baryon number, strangeness and charm) are equal to zero for this 
particle. Therefore, its production is not influenced by the conservation laws of quantum 
numbers. For suffiently high collision energies, the effect of the strict energy-momentum 
conservation in the statistical model formulation can be neglected. Consequently, the 
J/ip production can be calculated in the grand canonical approximation and, therefore, 
its multiplicity is proportional to the volume, V, of the matter at hadronization. Thus, 
the statistical yield of J/ip mesons at hadronization is given by 



where j = 1 and = 3097 MeV are the spin and the mass of the J/ip meson and Th 
is the hadronization temperature. The previously mentioned results of the analysis of 
hadron yield systematics in elementary and nuclear collisions within the statistical ap- 
proach indicate that the hadronization temperature Th is the same for different colliding 
systems and collision energies. This reflects the universal feature of the hadronization 
process. 

The total entropy of the produced matter is proportional to its volume. As most of the 
entropy in the final state is carried by pions, the pion multiplicity is also expected to be 
proportional to the volume of the hadronizing matter. Thus the scaling property (|l]) fol- 
lows directly from the hypothesis of statistical production of J/ ip mesons at hadronization 
and the universality of the parameter Th- 

Since elements of hadronizing matter move in the overall center of mass system the vol- 
ume V in Eq. (0) characterizes in fact the sum of the proper volumes of all elements in 
the collision event. 

The hypothesis of statistical production of J/ijj mesons at a constant hadronization 
temperature Th leads to the prediction of a second scaling property of the J/ ip multiplicity. 
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namely: 

^ consti^s) (3) 

[n ) 

which should be valid for sufficiently large cm. energies, a/s. This scaling property is 
illustrated in Fig. 2 which shows the ratio {J/ip) / {h~) as a function of ^/s for proton- 
nucleon interactions. The experimental data on J/ip yields are taken from a compilation 



given in [Q. The values of {h~) are calculated using a parameterization of the experi- 
mental results as proposed in |]T^. 

Onwards from the CERN SPS energies, ^/s = 20 GeV, the ratio {J/ilj)/{h^) is approx- 
imately constant, in line with the expected scaling behavior (^. The rapid increase of 
the ratio with collision energy observed below ^/s = 20 GeV should be attributed to a 
significantly larger energy threshold for the J/ ip production than for the pion production. 
In terms of the statistical approach the effect of strict energy-momentum conservation 
has to be taken into account by use of the microcanonical formulation of the model. 

The statistical J/ip multiplicity (^) depends on two parameters, and V. In general 
the calculation of the hadron yields in the statistical model should take into account the 
conservation of charges and resonance feeddown contributions. However, a simple way 
to estimate of the crucial temperature parameter in Eq. (|^) from the experimental data 
is possible, provided that we find a second hadron which has the properties of the J/ip 
meson i.e. it is neutral, unfiavored and stable with respect to strong decays. The best 
candidate is the t] meson. The multiplicity of t] mesons seems to obey also the scaling 
properties (0) and (^. We note, however, that the data on rj production are scarce. The 
independence of the (r]) / (tt^) ratio on the collision energy was observed quite a long time 
ago ||2^ . Recent data on t] production in central Pb+Pb collisions at the CERN SPS 



suggest that the (77) / (vr^) ratio is also independent of the size of the colliding objects. In 
order to illustrate this scaling Fig. 3 shows the (t]) / (tt") ratio as a function of the number 



of interacting nucleons for inelastic p+p |22| and S+S [23| interactions and for central 



Pb+Pb |T[ collisions at the CERN SPS energies (158-400 A GeV). 

From the ratios, (J/ip)/ (h^) and {r])/{'K^), presented in Fig. 1 and Fig. 3 we esti- 
mate a mean ratio {J/'ip)/{i]) = (1.3 ± 0.3) ■ 10~^. Here we use the experimental ratio 
(tt^)/ {h~) = 1 in N+N interactions [^. Under the hypothesis of the statistical produc- 
tion of J/tp and rj mesons at hadronization the measured ratio can be compared to the 
ratio calculated using Eq. (0): 

{J/^) ^ o f^A^^^ cxpf ""^"""^ ) (4) 



where = 547 MeV is the mass of the 77 meson. This leads to an estimate of the 
hadronization temperature, Th = 170 ± 2 MeV. Assuming a maximum 50% uncertainty 
on the {J/ip)/{ri) ratio due to the contribution from resonance decays^ we obtain an 
estimate of an additional systematic error on Th of about 7 MeV. A graphical solution of 

^ An estimate of the fraction of the 77 yield from decays of heavy hadrons in p+p interactions is about 
50% ||l^, which is close to the measured fraction of J/i/j yield originating from decays (30-50%) ||2^, |^. 
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Eq. is shown in Fig. 4 which illustrates the high sensitivity of the estimate of the Th 
parameter by using the {J /ip) / ijj) ratio. This is due to the large difference between mass 
of the J/ip and the r] mesons. 

In summary, we show that the J /if) production in hadronic and nuclear collisions can 
be understood assuming that a dominant fraction of J /tl) mesons is produced at hadroniza- 
tion according to the available hadronic phase space. The estimate of the hadronization 
temperature based on J/ip multiplicity, Th ~ 170 MeV, agrees well with the values of 
the temperature parameter obtained from the analysis of the hadron yield systematics in 
e"*" + e~, p + p, p + p interactions and nucleus-nucleus collisions. 

If the new interpretation of the J/ip data presented in this letter is correct, it may 
have several important implications: 

1. J/ip yields are not sensitive to the state of high density matter created at the early 
stage of A+A collisions because their production takes place at hadronization. 

2. The creation of the J/ip mesons is due to the straight thermal production at 
hadronization and not due to the coalescence of cc quarks produced before hadronization. 
Therefore the yield of J / ip mesons is independent of the production of open charm, which is 
carried mainly by the D mesons in the final state. The D meson multiplicity is determined 
by the number of cc quark pairs created in the early parton stage before the hadronization. 

3. Due to the large mass of J/ip mesons the data on their production are very sensitive 
to the value of the hadronization temperature and therefore allow for a precise study of 
the hadronization process. 
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Figure 1: The ratio of the mean muhipUcities of J/ip mesons and negatively charged 
hadrons for inelastic nucleon-nucleon (square) and inelastic 0+Cu, 0+U, S+U and 
Pb+Pb (circles) interactions at 158 A-GeV plotted as a function of the mean number 
of participant nucleons. For clarity the N+N point is shifted from (Np) = 2 to (Np) — 5. 
The dashed line indicates the mean value of the ratio. 
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Figure 2: The ratio of the mean multiphcities of J/ if) mesons and negatively charged 
hadrons for inelastic proton-nucleon interactions as a function of the collision energy in 
the center of mass system. 
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Figure 3: The ratio of the mean multiphcities of rj mesons and 7r° mesons for inelastic 
p+p and S+S interactions and central Pb+Pb collisions at the CERN SPS energies as 
a function of the number of interacting nucleons. The result for Pb+Pb collisions was 
obtained in the rapidity interval 2.3-2.9. 
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Figure 4: The {J/ip)/{ri) ratio calculated under hypothesis of the statistical production 
of J/ip and 7] mesons at hadronization (solid line) as a function of the hadronization 
temperature. Band shown by dashed lines is drawn at ±a around the mean experimental 
value of the {J/ip)/ {rf) ratio. The dotted line indicates Th = 170 MeV. 
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